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Resumen 


Inicialmente se expone una clasificación para los modelos probabilísticos 
utilizados en el análisis de frecuencia de crecientes local. Enseguida se presenta un 
resumen de los aspectos teóricos y estadísticos asociados con el llamado modelo 
de analogía de difusión lineal, el cual constituye una distribución probabilística 
físicamente basada. Posteriormente se describe su aplicación numérica en cinco 
registros de gastos máximos anuales de la cuenca baja del río Pánuco. Por último, 
se analizan los resultados y se formulan las conclusiones, las cuales destacan la 
conveniencia de este modelo en estaciones hidrométricas de gran área drenada en 


ríos de planicie. 
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Introducción 


El análisis de frecuencia de crecientes (AFC) 
ha sido un tema de gran interés práctico y de 
mucha investigación dentro de la hidrología 
y el desarrollo de los recursos hidráulicos, 
principalmente debido a que permite el dimen- 
sionamiento de la infraestructura hidráulica, 
desde los vertedores de los embalses para 
brindar seguridad en éstos, hasta las obras de 
cruce (alcantarillas y puentes) y de protección 
(diques, 
Como resultado de tales tendencias, en el AFC se 


encauzamientos y rectificaciones). 
utilizan diversos enfoques y variados modelos 
matemáticos (Campos, 1999), los cuales, de 
manera general, se clasifican en locales o del 
sitio, y regionales. Por su parte, los modelos 
probabilísticos que se utilizan en el AFC local 
se pueden clasificar en los siguientes tres tipos 
(Strupczewski et al., 2001): (1) empíricos, (2) 
fenomenológicos y (3) físicamente basados. 

El enfoque empírico es actualmente el 
más común. Consiste en utilizar modelos o 
distribuciones de probabilidad que han sido 
propuestas para realizar el AFC, mismas que 
han demostrado su conveniencia o bondad de 
ajuste. Stedinger et al. (1993), y Rao y Hamed 


(2000) exponen con detalle esta aproximación. 
En el enfoque fenomenológico se utilizan 
axiomas probabilísticos, los cuales conducen a la 
formulación de un modelo matemático de una o 
más características de las crecientes. Son ejemplos 
de esta aproximación la modelación estocástica 
de crecientes y la aplicación de la teoría de la 
entropía. Este tratamiento no llegó a ser popular 
debido a su sofisticación matemática. Por 
último, en el enfoque de los modelos físicamente 
basados se emplean principios dinámicos de 
la generación de las crecientes para obtener 
un modelo probabilístico. En este contexto, 
Strupczewski et al. (2001) utilizaron los modelos 
de tránsito de gasto en canales, en los cuales no 
se hace una consideración explícita a los procesos 
hidrológicos que ocurren en los suelos de la 
cuenca para desarrollar una función de densidad 
de probabilidad con dos parámetros de ajuste. 

El objetivo de este trabajo consiste en 
exponer un resumen de la teoría del modelo de 
analogía de difusión lineal (MADIL) desde su 
concepción hasta las ecuaciones de su aplicación 
como distribución probabilística, incluyendo 
su aplicación y contraste en cinco estaciones 
hidrométricas de la cuenca del bajo río Pánuco o 
Región Hidrológica Número 26, de México. 
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Desarrollo 
Resumen de teoría 


Aceptando que el volumen unitario de la 
función impulso 5 de Dirac consiste de un 
número infinito de partículas o gotas, entonces 
la integral de la respuesta impulso, es decir: 


Jrs, £)- at (1) 


determina la probabilidad de que una partícula 
individual pase por la salida s durante el tiempo 
(0,7), donde h(s,t) es la función a la respuesta 
impulso en el tiempo t en la posición s. Además 
de esta interpretación probabilística, la función 
citada debe cumplir con los requerimientos 
de los modelos probabilísticos empleados en 
el AFC, como son: límite inferior no negativo, 
sesgo positivo y ser unimodal. 

Por otra parte, la predicción de las 
características de una creciente (gasto pico y 
volumen) en una sección aguas abajo de un 
río, con base en las condiciones conocidas del 
flujo o hidrograma en una sección de aguas 
arriba, es el problema fundamental del tránsito 
de crecientes. Usando una linealización de las 
ecuaciones de Saint Venant (Dooge y Harley, 
1967; Napiórkowski, 1992) se obtiene un 
modelo lineal físicamente basado con cuatro 
parámetros dependientes de las características 
hidráulicas del tramo de cauce, cuya solución 
aún es compleja. Se desarrollaron dos versiones 
simplificadas: el modelo de analogía de difusión 
lineal o modelo convectivo-difusivo lineal, 
y el modelo de flujo rápido lineal, los cuales 
corresponden a condiciones límites del flujo 
establecidas con el número de Froude de cero y 
de uno, respectivamente. La respuesta impulso 
del primer modelo está dada por (Strupczewski 
et al., 2001): 


pel 
Hs h= art) | (2) 


E exp 
V4TD 1? 4D-t 


donde s es la longitud de tramo de canal, t es 
el tiempo, u es la velocidad convectiva y D es 
la difusividad hidráulica. Tanto u como D son 
funciones de las características del canal y 
del flujo en la condición de referencia de flujo 
estable. Designando a s/V4D=0, a s-u/4D=P 
y at como x, se obtiene la función de densidad 
de probabilidad con dos parámetros de ajuste 
del modelo de analogía de difusión lineal 
(MADIL): 


Ax:0,P) == 


a. es el parámetro de escala con dimensión de 
raíz cuadrada de las unidades de la variable 
aleatoria x y P es el parámetro de forma 
adimensional. Este modelo probabilístico ha 
sido estudiado desde finales de la década de 
los años cincuenta por diversos matemáticos y 
estadísticos, siendo denominado distribución 
Gaussiana inversa. 


Principales propiedades estadísticas 


Los coeficientes de variación, asimetría y 
kurtosis son (Strupczewski et al., 2001): 


1 
3 EL 4 
Cv 2B (4) 
cate (5) 
2B 
5 
Cha zp =3(5Cv? +1) (6) 
la moda es: 
Mo=% E 1 (7) 
ABRI WUS 


Cuando Cv>0, la distribución tiende a 
ser simétrica, como los modelos Log-Normal 
y Gamma. La ecuación (5) establece que en el 
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MADIL, el cociente Cs/Cv es igual a tres; en 
cambio, en la distribución Gamma es de dos, y 
en la Log-Normal existe esta relación: Cs = 3Cv 
+ Cv?, Estas tres distribuciones y la Weibull han 
sido contrastadas por Mitosek et al. (2006) en 
varios ríos de Polonia, con registros de gastos 
máximos anuales del orden de setenta años, 
encontrado como mejores opciones el MADIL 
y la distribución Log-Normal. Para todas las 
distribuciones citadas se han desarrollado 
que incluyen dos 
parámetros de forma, con la idea de buscar 


versiones matemáticas 


un mejor ajuste en su porción o cola derecha 
(Strupczewski et al., 2008). 


Estimación de parámetros y de su dispersión 
Los métodos de momentos (MMO) y de 


(MMV) conducen 
respectivamente a las expresiones siguientes 


máxima verosimilitud 


para los parámetros de escala y forma 
(Strupczewski et al., 2001): 


BR E pel 

a Ñ 
z 

is a S Ez 1) da 


en donde x es la media aritmética y H, la media 
de los recíprocos; es decir: 


n 

¿Y x, (12) 
i=1 
n 

H - > (13) 


siendo x, los gastos máximos anuales (m*/s) de 
la serie o muestra que se procesa, cuyo número 
de datos es n. 


Las ecuaciones (10) y (11) del MMV están 
basadas en las medias aritmética y de recíprocos, 
por lo tanto sus estimaciones no son sensibles a 
los grandes valores de x en la muestra (outliers). 
Lo anterior significa que este método conduce 
a un buen ajuste en la porción izquierda de la 
distribución; en cambio, con el MMO se logra un 
mejor ajuste en parte extrema de los datos y por 
ello es el preferido (Strupczewski et al., 2001). 

Las dispersiones de los parámetros de ajuste 
D(a) y D(B), expresadas como error estándar, 
son en el MMO las siguientes: 


D*(o)= E (acots 1) (14) 
pp (7C0?+2) (15) 


Se observa que en ambos casos tales 
dispersiones dependen directamente de la 
variabilidad de los datos e inversamente de la 
longitud de la muestra. 


Estimación de predicciones 
El valor de x asociado con una determinada 
probabilidad de excedencia p (x,) o predicción 
se obtiene integrando la ecuación (3), lo cual 
conduce a: 
(16) 
en la cual q, es el límite superior de la siguiente 
integral: 

Op 


A] E se 


La solución de tal integral se obtuvo 


(17) 


mediante la cuadratura de Gauss, utilizando 
el polinomio de Legendre aplicado a cualquier 
intervalo de integración y mediante una 
fórmula de diez puntos (Nakamura, 1992; 
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Campos, 2003). Algunos valores de q, para las 
magnitudes de f y de p que ocurren y se usan 
en la práctica se exponen en el cuadro 1. 
Strupczewski et al. (2001) presentan las 
ecuaciones para estimar la dispersión en 
las predicciones Dx), así como sus valores 


relativos [vz -D(xp) / ña] en forma tabular para 


las magnitudes de B y de p que ocurren y se 
usan en la práctica. "lales valores relativos se 
reproducen en el cuadro 2. 


Aplicaciones numéricas 


Strupczewski et al. (2001) encontraron que el 
ajuste del MADIL fue notablemente mejor en 


los registros de las estaciones hidrométricas con 
grandes áreas de cuenca (> 10 000 km?), es decir, 
en ríos de planicie, donde el modelo convectivo- 
difusivo lineal es una buena aproximación 
para modelar el tránsito de crecientes. Por ello, 
las aplicaciones numéricas del MADIL que 
serán expuestas con propósitos descriptivos se 
realizaron en las estaciones de aforos finales de 
la Región Hidrológica Número 26 (Pánuco), las 
cuales se ubican en la planicie costera del golfo 
de México. 

En el cuadro 3 se exponen las caracterís- 
ticas generales de las cinco estaciones 
hidrométricas procesadas y las propiedades 
estadísticas de sus registros de gasto máximo 


Cuadro 1. Magnitudes del límite q, (adimensional) para los valores de la probabilidad 


de excedencia (p) y de B indicados. 


Te (años) | p (9%) Valores del parámetro de forma $ 

0.50 1.00 150 2.00 2.50 3.00 3.50 4.00 4.50 5.00 
2 50 0.8602 | 1.1150 1.3209 1.4986 | 1.6573 020) [155 240510 | 21700 | 22505 
O 20 0.5877 | 0.8414 1.0457 | 1.2219 1.3793 1.5230 1.6559 1.7804 1.8977 | 2.0090 
10 10 0.4830 | 0.7282 | 0.9272 | 1.0999 1.2546 1.3961 SS 1.6503 1.7664 | 1.8766 
25 4 0.3960 | 0.6276 | 0.8185 | 0.9856 | 1.1362 1.2744 | 1.4030 | 1.5238 1.6380 | 1.7466 
50 2 OSSI0 1015725 0.7571 0.9199 1.0672 1.2029 148295 1.4484 1.5611 1.6684 
100 1 0.3168 | 0.5284 | 0.7073 | 0.8658 1.0099 1.1430 | 1.2674 | 1.3847 1.4958 1.6018 
500 0.20 | 0.2617 | 0.4539 0.6201 0.7694 | 0.9064 MOSS MELOSS 1.2667 1.3744 | 1.4774 
1 000 0,10 | 0.2446 | 0.4295 0.5909 | 0.7366 | 0.8707 | 0.9958 1.1135 1.2250 1.3313 1.4331 
5 000 0.02 | 0.2142 | 0.3845 0.5357 | 0.6737 | 0.8017 | 0.9217 | 1.0352 | 1.1431 1.2461 1.3450 
10 000 0.01 | 0.2039 | 0.3688 0.5160 | 0.6511 0.7766 | 0.8946 | 1.0063 1.1127 1.2144 | 1.3122 

* Periodo de retorno. 


Cuadro 2. Errores relativos Jn -D (»,) /xp (adimensionales) para los valores de la probabilidad 


de excedencia (p) y de P indicados. 


T+* (años) | p (9%) Valores del parámetro de forma f 
0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 

2 50 1.1271 | 0.7147 [0.5692 | 0.4897 | 0.4373 | 0.3994 | 0.3701 | 0.3466 | 0.3272 | 0.3107 
a 20 1.0013 | 0.7644 | 0.6432 | 0.5654 | 0.5101 | 0.4682 | 0.4351 | 0.4079 | 0.3853 | 0.3661 
10 10 1.2634 | 0.9470 | 0.7851 | 0.6837 | 0.6128 | 0.5598 | 0.5183 | 0.4847 | 0.4568 | 0.4332 
20 5 15450 | 1.1238 | 0.9206 | 0.7963 | 0.7106 | 0.6472 | 0.5979 | 0.5583 | 0.5253 | 0.4976 
50 2 1.8572 | 1.3224 | 1.0748 | 0.9660 | 0.8242 | 0.7495 | 0.6917 | 0.6451 | 0.6067 | 0.5744 
100 1 2.0468 | 1.4464 | 1.1731 | 1.0099 | 0.8989 | 0.8172 | 0.7540 | 0.7032 | 0.6614 | 0.6262 
500 0.20 2.3677 | 1.6684 | 1.3552 | 1.1685 | 1.0413 | 0.9479 | 0.8757 | 0.8177 | 0.7697 | 0.7291 
1000 0.10 2.4703 | 1.7424 | 14178 | 1.2242 | 1.0925 | 0.9954 | 0.9202 | 0.8597 | 0.8097 | 0.7674 


* Periodo de retorno. 


Campos-Aranda, D.F., Descripción de un modelo físicamente basado para el análisis de frecuencia de crecientes 


anual. Las series de las estaciones "Tamesí y 
Pánuco han sido expuestas por Campos (2007) 
y la de Magiscatzin II se presentó en Campos 
(2006). Las muestras de las estaciones El Pujal 
y Tempoal proceden del sistema BANDAS 
(IMTA, 2002). Todos los registros citados se 
tienen en el cuadro 4. 

Por otra parte, en el cuadro 5 se exponen los 
resultados de las ecuaciones (8), (9), (14) y (15); 
es decir, del cálculo de los parámetros de ajuste 
y sus dispersiones. 


Contraste y error estándar de ajuste 


A los cinco registros hidrométricos procesados 
(cuadro 4) se les aplicaron las distribuciones 
Log-Pearson tipo III (LP3), a través del método 
de momentos en sus dominios logarítmico y 
real (Bobée y Ashkar, 1991), así como la General 
de Valores Extremos (GVE), por medio de 
los métodos de momentos, sextiles, máxima 
verosimilitud y momentos L (Campos, 2001), 
seleccionando el ajuste que condujo al menor 
error estándar de ajuste (EEA), definido como 
(Kite, 1977): 


(18) 


en donde n es el número de datos del registro; 
x, los gastos máximos anuales observados, 


ordenados de menor a mayor; X, los gastos 


máximos estimados con el modelo probabilístico 
que se prueba, para una probabilidad de no 
excedencia, definida con la fórmula Weibull 
(Benson, 1962): 


P(x<x)= 


n+1 (19) 


en la cual X es la variable aleatoria, cuyo límite es 
x¿ m es el número de orden del dato (x,), con uno 
para el menor y 1 para el mayor. Finalmente, np 
es el número de parámetros de ajuste: np=2 para 
el MADIL y np=3 para los modelos LP3 y GVE. 

Como el MADIL (ecuaciones (16) y (17) no 
tiene solución inversa; es decir, que no se puede 
estimar directamente el valor de x asociado con 
una determinada probabilidad de no excedencia 
F(x), entonces, para el cálculo del EEA, se utilizó 
el procedimiento de bisección sugerido por 
Campos (1989); para ello se proponen valores 
iniciales inferior (xi=0.0001:x) y superior 
(xs=10-x,); se obtiene su media aritmética (%;); 
con tal magnitud se obtiene q, (ecuación (16)) y 
después F(%)=p con la ecuación (17); se compara 
este último valor con el de la ecuación (19), 
relativo al x, y si su diferencia absoluta no es 
menor que 0.00001, que es la tolerancia aceptada, 
se sigue el método de bisección. 


Análisis de los resultados 


En el cuadro 6 se han concentrado las 
predicciones obtenidas con los tres modelos 
probabilísticos aplicados, en las cuales se 
observa gran similitud en los periodos de 


Cuadro 3. Características generales de las estaciones hidrométricas y de sus registros 


de gastos máximos anuales utilizados. 


S Área A Núm. 
Nombre Clave Corriente Registro y Cv Cs Cs/Cv 
(km?) años 
Tempoal 26248 Río Tempoal 5275 1954-2002 49 0.7432 1.3988 1.882 
Magiscatzin II 26249 Río Guayalejo 10 968 1954-2002 49 0.7519 bso 2.013 
Tamesí 26416 Río Tamesí 14 923 1973-2002 30 0.4687 1.0725 2.288 
El Pujal 26272 Río Tampaón 23373 1954-2002 49 0.6376 1.4190 2.226 
Pánuco 26424 Río Pánuco 65 577 1972-2002 31 0.4119 1.5790 3.833 


Simbología: Cv: coeficiente de variación; Cs: coeficiente de asimetría. 
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Cuadro 4. Gastos máximos anuales (m?/s) disponibles en las estaciones hidrométricas indicadas 
de la cuenca del Bajo Río Pánuco. 


Año Tempoal Magiscatzin II Tamesí El Pujal Pánuco 

1954 2 110.0 524.0 3 100.0 

1955 6 000.0 4. 500.0 6 052.0 

1956 4 424.0 1 150.6 2 395.0 

1957 449.0 261.6 564.3 

1958 4 100.0 491.3 3 064.0 

1959 1 507.6 1 384.0 685.0 

1960 1277.0 509.9 854.0 

1961 852.9 112399 1 449.3 

1962 739.2 1 589.0 484.4 

1963 1 800.0 363.6 1516.6 

1964 748.0 513.8 428.0 

1965 EDAD 851.5 1 125.8 

1966 1778.0 3737.5 2 230.0 

1967 2 245.0 2 124.2 2 860.0 

1968 1 145.0 908.1 983.0 

1969 1 948.0 633.0 2 845.3 

1970 1 418.0 3 403.6 2 630.0 

1971 1 630.0 2311.6 1 467.6 

1972 989.0 2 286.0 1 585.0 2 880.0 

1973 1 668.0 2 997.4 2 700.0 2 064.0 3 234.0 

1974 4 950.0 796.0 1710.0 4 140.0 7 300.0 

1975 4 040.0 1 460.3 1 102.0 3 038.0 4 138.0 

1976 1275.0 5 510.0 641.0 4 902.0 3 886.0 

1977 514.0 2 600.0 2 148.0 2 683.0 1 968.0 

1978 3725.0 967.0 1 157.0 3 056.0 3 471.0 

1979 1 655.9 778.4 1 147.0 1 301.0 3 525.0 

1980 1162.0 429.5 449.5 1 142.9 2 753.0 

1981 2 020.0 934.5 1 170.0 1 583.3 3 240.0 

1982 539.6 349.5 427.5 364.6 829.9 

1983 868.0 996.0 1 180.0 171.5 2 631.0 

1984 4 030.0 584.4 1 124.0 1987.1 4 236.0 

1985 1 882.0 1071.1 1 140.0 1 346.0 2 906.0 

1986 476.0 742.0 674.0 946.0 1 641.0 

1987 1765.0 2 040.0 1 525.0 896.0 2775.00 

1988 3 265.0 1 190.0 1256.0 1 625.0 3 170.0 

1989 649.0 653.0 925.0 696.0 1 530.0 

1990 1 611.0 534.0 850.0 2 760.0 3 445.0 

1991 3 532.0 2 500.0 1 750.0 2 270.0 3 619.0 

1992 2 291.0 3111.9 1918.1 2 165.0 2 890.0 

1993 6 120.0 2 231.6 2 387.0 4 760.5 5 400.0 

1994 1 133.0 1 648.2 1 202.6 1 184.6 2 325.0 
e 1995 742.0 TAS 965.9 1 404.9 1976.0 
5 1996 683.0 1 401.5 1 063.4 2 097.2 2 546.0 
% 1997 905.0 1 036.8 990.7 ISSiS 2 334.0 
El 1998 1266.9 794.5 870.6 1 249.8 2 486.6 
E 1999 2 693.7 1921.6 1 092.9 1 244.1 2 701.4 
S 2000 641.2 2 931.6 2 577.7 1 128.1 1 818.0 
E 2001 1 847.9 1 059.2 917.9 1 296.2 2 306.4 
$ 2002 926.4 991.5 689.5 976.7 1 899.4 
3 
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Cuadro 5. Parámetros de ajuste y su dispersión en el MADIL ajustado por momentos en las cinco 


estaciones hidrométricas indicadas de la Región Hidrológica Número 26. 


QL D(o)/o. DB)! 
Nombre (m/s) Día) % D(B) es B 
Tempoal 41.8569 0.9053 6.8920 16.47 0.3132 34.60 
Magiscatzin II 36.6847 0.8844 6.0846 16.59 0.3084 34.87 
Tamesí 53.5143 2.2757 8.8989 16.63 0.7815 34.34 
El Pujal 48.4409 1.2301 7.2899 15.05 0.3868 31.45 
Pánuco 93.4453 2.9468 14.5784 15.60 0.9449 32.07 
Cuadro 6. Predicciones de gasto máximo (m? /s) en las cinco estaciones hidrométricas indicadas 
de la Región Hidrológica Número 26 (Pánuco). 
Periodos de retorno en años 
Nombre FDP? EEA! (nY/s) 
2 5 10 25 50 100 500 1000 5000 | 10000 
Tempoal 
Log-Pearson tipo 
m 285 1478 | 2700 | 3770 | 5543 | 6972 | 8739 | 14005 | 16874 | 25301 | 29834 
GVE? 303 (15209 4216979 ESI0704 Mos 18 [MOJ5558 ESEL69A MESAS 7 SOSA 226398 291561 
MADIL* 285 1526 | 2749 | 3721 | 5074 | 6146 | 7254 9938 | 11133 | 13981 | 15234 
Magiscatzin II 
Log-Pearson tipo 
m 128 1170 | 2158 | 299 | 4274 | 5391 | 6655 | 10245 | 12115 | 17350 | 20041 
GVE? 163 1205 | 2144 | 2914 | 4094 | 5150 | 6379 | 10089 | 12154 | 18409 | 21890 
MADIL* 177 119 | 2165 | 2939 | 4020 | 4878 | 5765 7918 8877 | 11165 | 12171 
Tamesí 
Log-Pearson tipo 
m 135 1143 | 1673 | 2040 | 2516 | 2878 | 3247 4 136 4.536 5512 | 5953 
GVE? 135 1125 | 1640 | 2022 | 2557 | 2995 | 3468 4725 5 345 6999 | 7816 
MADIL* 135 1136 | 1667 | 2033 | 2500 | 2849 | 3195 3999 4 346 5155 | 5505 
El Pujal 
Log-Pearson tipo 
m 189 1613 | 2700 | 3501 | 4581 | 5427 | 6301 8 452 9434 | 11847 | 12944 
GVE? 190 1619 | 2656 | 3445 | 4575 | 5522 | 6564 9428 | 10885 | 14902 | 16 947 
MADIL* 203 159 | 2657 | 3457 | 4538 | 5375 | 6228 8 264 9162 | 11284 | 12213 
Pánuco 
Log-Pearson tipo 
m 339 2734 | 3812 | 4554 | 5519 | 6257 | 7010 8 845 9 681 11744 | 12690 
GVE? 349 2757 | 379% | 4509 | 5442 | 6158 | 6888 8 658 9 455 so [12255 
MADIL* 343 2734 | 3839 | 4575 | 5499 | 6178 | 6848 8 383 9040 | 10561 | 11216 
Simbología: 


* Error Estándar de Ajuste. 
Función de Distribución de Probabilidades. 
3General de Valores Extremos. 

“Modelo de Anaogía de Difusión Lineal. 
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retorno bajos (<50 años); en cambio, en los 
intervalos de recurrencia extrema (>500 años) 
se tienen diferencias significativas en Tempoal 
y Magiscatzin Il que son las estaciones 
hidrométricas con los cocientes Cs/Cv más 
alejados de tres. 

Al respecto, Strupczewski et al. (2001) 
encontraron en 39 ríos polacos un cociente 
Cs/Cv medio de 2.52; esto es, un valor más 
cercano al MADIL que al de las distribuciones 
Gamma y Log-Normal; por ello, tal modelo fue 
el más conveniente en 27 de tales registros. 

En el resto de las estaciones, las predic- 
ciones del MADIL coinciden con las de las 
distribuciones LP3 y GVE, que son modelos 
más flexibles al tener tres parámetros de 
ajuste, mismos que han demostrado su bondad 


de ajuste y por ello su aplicación ha sido 


establecida como precepto en Estados Unidos 
e Inglaterra, respectivamente. 

Respecto al error estándar de ajuste (EEA), 
en general son del mismo orden de magnitud 
en los tres modelos probabilísticos aplicados. 
El MADIL en la estación El Pujal conduce al 
mayor EEA, sin embargo sus predicciones son 
coincidentes con las de la distribución LP3. 

En la figura 1 se muestra el ajuste del 
MADIL a los datos de la estación hidrométrica 
Pánuco; además, se indica la dispersión de las 
predicciones, estimadas para B=3 con base en 
los valores del cuadro 2. 


Conclusiones 


En estaciones hidrométricas con gran área 
drenada ubicadas en ríos de planicie, en los 
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Figura 1. Contraste gráfico entre el modelo de analogía de difusión lineal y los gastos máximos 


observados en la estación hidrométrica Pánuco. 
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cuales la aproximación lineal para el tránsito de 
crecientes es válida y donde además el cociente 
entre los coeficientes de asimetría y de variación 
de su registro de gastos máximos anuales sea 
cercano a tres, es muy probable que el MADIL 
sea la distribución de probabilidades más 
conveniente. 

Este modelo probabilístico físicamente 
basado tiene la ventaja de ser fácilmente aplica- 
ble y cuya evaluación de sus parámetros de ajuste 
utiliza dos propiedades estadísticas básicas: 
la media aritmética y la desviación estándar a 
través del Cv, esto implica una estimación más 
confiable. 
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Abstract 


CAMPOS-ARANDA, D.F. Description of a physically based model for flood frequency 
analysis. Water Technology and Sciences, formerly Hydraulic engineering en Mexico 
(in Spanish). Vol. [, no. 2, April-June, 2010, pp. 113-122. 


First, a classification is presented for the probabilistic models used in local flood frequency 
analysis. Then, a summary is given of the theoretical and statistical aspects associated with the 
so-called linear difusion analogy model, which is a physically based probabilistic distribution. 
Next, its numerical application in described in five records of annual maximum flows of the 
lower Pánuco River watershed. Finally, results are analyzed and conclusions drawn, which 
point to the suitability of this model in hydrometric stations with a large drainage area in 
lowland rivers. 


Keywords: linear model, diffusion model, probability distributions, method of moments, 
standard error of fit, Pánuco River. 
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